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INTRODUCTION: PROGRAM OBJECTIVES 


The goal of this program, supported by NASA-Lewis Research Center from October 25, 1985 
to February 28, 1989, was elucidation of the chemical reaction mechanism responsible for soot 
formation in hydrocarbon flames and development of computer models for practical applications. 
This program is a continuation of the work, initiated under the previous NASA Grant No. NAG 3- 
477, which resulted in discovery of the H-abstraction/C 2 H 2 -addition mechanism responsible for 
the formation and growth of polycyclic aromatic hydrocarbons (PAHs). 

The primary focus of the work reported here was on kinetic modeling. The specific objectives 
were: detailed modeling of soot formation in premixed flames, elucidation of the effects of fuel 
structure on the pathway to soot, and development of a numerical technique for accurate modeling 
of soot particle coagulation and surface growth. We successfully completed all the above tasks, 
which are briefly summarized below. 

DETAILED MODELING OF SOOT FORMATION IN FLAMES 

A detailed modeling study of the formation of PAHs in a burner-stabilized low-pressure 
sooting 23.6 % acetylene - 21.4 % oxygen - argon laminar premixed flame of Bockhorn and co- 
workers 1 was performed. The model predicted the correct orders of magnitude and relative 
appearances of the concentration peaks, but overstated the decline of the species concentrations in 
the post-flame zone. Imprecise knowledge of the thermochemical data and unknown details of the 
oxidation of hydrocarbon radicals are the reasons identified for the latter. The main reaction 
pathways for cyclization and growth of PAHs and the results of the sensitivity tests are in close 
agreement with those of the previous modeling study of acetylene oxidation under shock-tube 
conditions. It was identified that an additional factor important in the flame environment is the 
diffusion of hydrogen atoms from the main reaction zone into a cooler preheat region. 

ELUCIDATION OF REACTION PATHWAYS TO SOOT 
Effect of Fuel Structure 

The influence of fuel structure on chemical reaction pathways leading to soot formation was 
investigated using detailed kinetic models for the conditions of shock-tube pyrolysis of 1,3- 
butadiene, benzene and ethylene. It was found that reaction pathways to soot relax to the 
previously identified H-abstraction/C 2 H 2 *addition mechanism^ in each case. Fuel molecular 
structure influences the growth process only at its early stages, first, by providing more efficient 
reaction partners for the formation of aromatic species and, second, by affecting generation of 
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hydrogen atoms. The model explains the relative sooting tendencies of these fuels for the shock- 
tube pyrolysis conditions 3 ' 5 and agreed with recent observations in laminar diffusion flames. 6 ' 8 

Experimental and Modeling Study of Soot F orma tion in Mixtures sil 

Hydrocarbons 

Soot formation in hydrogen-, allene-, vinylacetylene-acetylene and acetylene-, allene-, 
vinylacetylene- and 1,3-butadiene-benzene argon-diluted mixtures was studied behind reflected 
shock waves by monitoring the attenuation of a 632.8 nm He-Ne laser beam. The experiments 
were performed at temperatures 1500-2490 K, pressures 1.3-3. 1 bar, and total carbon atom 
concentrations (2.0-6.8) x 10 17 atoms/cm 3 . The experimental observations were interpreted using 
detailed chemical kinetic modeling. The reaction mechanism used in the simulations was 
composed of 520 reversible reactions and 180 species. 

The addition of hydrogen to acetylene resulted in the decrease of soot production. Pronounced 
synergistic effects were observed for soot amounts in all binary hydrocarbon mixtures studied. 
However, the same experimental results expressed in terms of soot yields showed much smaller 
effects and in some cases, like in benzene-acetylene mixtures, the negative effects were obtained. 
These observations indicate that a) competition between promoting and suppressing factors takes 
place when individual fuels are mixed together, and b) one must distinguish between various 
experimental characteristics in the analysis of sooting trends in fuel mixtures. 

The computational results, which were found to be in qualitative agreement with experiment, 
revealed that the main factor affecting PAH (and soot) formation in the hydrocarbon mixtures is the 
acceleration of acetylene-addition reactions. In mixtures of aliphatic fuels, e.g. butadiene- 
acetylene, most of these reactions are those of the dominant ring formation and growth pathway, 
while in benzene-additive mixtures they play different, competing roles. The main route to PAH in 
benzene pyrolysis is the formation of biphenyl by addition of phenyl to benzene followed by 
sequential addition of acetylene to form pyrene. The increase in the rate of acetylene-addition 
reactions of this route enhances the production of soot. At the same time, reaction of acetylene 
with phenyl forming phenylacetylene removes phenyl radicals from a more efficient ring-forming 
pathway, addition to benzene molecules, which counteracts the promoting factor. In the mixtures 
of aliphatics the route to soot is that via phenylacetylene and thus no competition for phenyl exists. 

The influence of vinylacetylene and butadiene on PAH formation from benzene appeared to be 
mainly in the production of acetylene. The results obtained in benzene-allene mixtures were 
interpreted in an analogous way. It was suggested that the first-ring cyclization in allene mixtures 
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proceeds not via a direct formation of benzene, as usually assumed, but by sequential addition of 
two C2H2 molecules to C3H3 to form benzyl. 

The results of this study, both experimental and computational, provide further evidence for the 
importance of interaction between aromatic rings and acetylenic species, and that it is acetylene 
whose reactions appear to be of particular importance. 

Soot Formation From Vinylacetylene 

Soot formation in vinylacetylene, and vinylacetylene-oxygen argon-diluted mixtures was 
studied behind reflected shock waves by monitoring the attenuation of a 632.8-nm He-Ne laser 
beam. The experiments were performed at temperatures of 1600-2500 K, pressures of 2.08-3.09 
bar, and total carbon atom concentrations of (1.99-2.05)xl0 17 atoms/cm 3 . The experimental 
results obtained in pyrolysis of vinylacetylene are similar to those of acetylene, both in the order of 
magnitude of the soot yield and the shape of its temperature dependence. It is important that there 
is a secondary increase in the amount of soot formed from vinylacetylene, after about 2300 K, 
similarly to the phenomenon observed in the pyrolysis of acetylene. The addition of oxygen to 
vinylacetylene shifts the “soot bell” to lower temperature and, distinct from all other hydrocarbons 
studied in our laboratory (acetylene, 1,3-butadiene, allene, benzene and toluene), accelerates the 
production of soot with reaction time. Thus, it is shown that the soot amounts are smaller in a 
vinylacetylene-oxygen mixture than in a comparable vinylacetylene mixture at a reaction time of 0.5 
ms, but with the increase in reaction time this ratio is reversed. 

A detailed chemical kinetic simulation of soot formation was performed for the conditions of 
vinylacetylene pyrolysis. The reaction model used in this study was the one developed for the 
analysis of soot formation in binary hydrocarbon mixture. The computational results explain the 
experimentally observed similarity in soot formation from vinylacetylene and acetylene. 
Vinylacetylene is hard to decompose. When it does, which takes place at high temperatures, it 
decomposes quickly to form acetylene, which then explains the noted similarity. Oxygen promotes 
the decomposition of vinylacetylene, as it does in the case of other fuel molecules. However, in 
the case of vinylacetylene there is a high concentration of C4H3 radicals. Having them and 
acetylene molecules in high concentrations simultaneously accelerates the production of aromatic 
species and soot, which explains the effect of oxygen on soot formation observed in this study. 
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Soot Formation in Benzene-Alcohol Mixtures 


Soot formation in benzene-methanol and benzene-ethanol argon-diluted mixtures was studied 
behind reflected shock waves by monitoring the attenuation of a He-Ne laser beam. The 
experiments were performed at temperatures 1580-2250 K, pressures 2.0-3.0 bar, and total carbon 
atom concentrations (2. 0-2. 7) x 10^ atoms/cmA The results obtained indicate that the addition of 
alcohol suppresses the formation of soot from benzene at all temperatures and the reduction in soot 
yields is increased with the amount of alcohol added. The analysis of the results indicates that the 
suppression effect is probably due to oxidation of soot and soot precursors by OH and removal of 
hydrogen atoms by alcohol and water molecules. 

Spheroidal Carbon Clusters in Soot Formation 

It has been recently proposed^ that primary soot particles, being produced rapidly in flame 
environment and reaching typically an average diameter of 20 to 40 nm, are formed by a spheroidal 
shell growth of partially closed carbon clusters. The analysis of such a possibility in light of our 
chemical kinetic model for PAH and soot formation along with the results of electron diffraction 
analysis obtained by Ebert at Exxon led us to conclude that the importance of the carbon clusters is 
very unlikely in soot formation. 

ANALYSIS OF THE DRIVING FORCE OF PAH FORMATION 

A theoretical analysis of the kinetic factors affecting the production rate of polycyclic aromatic 
hydrocarbons in high-temperature pyrolysis and combustion environments was developed. I he 
analysis is based on a lumped kinetic model, representing polymerization-type growth by one 
irreversible and two reversible steps, which underlies the H-abstraction/C 2 H 2 -addition reaction 
mechanism identified in previous detailed kinetic studies. The analytical solution obtained for the 
steady-state assumption indicates that PAH production is not determined by a single parameter; 
several variables influence the process to different degrees depending on experimental conditions. 
At high temperatures, PAH growth is controlled by the superequilibrium of hydrogen atoms. At 
low temperatures and low H 2 concentrations, the PAH growth rate is proportional to the rate of the 
H-abstraction of a hydrogen atom from aromatic molecules; at low temperatures and high H 2 
concentrations, it is controlled by the thermodynamics of the H-abstraction and the kinetics of 
acetylene addition to aromatic radicals. The presence of oxygen mainly affects the small-molecule 
reactions during the induction period. 
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AEROSOL DYNAMICS 


A new numerical method to the solution of an aerosol dynamics has been developed. It is 
accurate, extremely economical, applicable to any combination of aerosol processes, and does not 
require prior knowledge of the particle size distribution function. The method is based on the 
description of the aerosol evolution by moments of the particle size distribution function. The 
underlying idea of the method of moments is that the properties one considers in most practical 
applications are determined by average measures; the history of individual particles can be ignored, 
and what is lost in resolution can be compensated by a significant increase in computational speed 
and a dramatic reduction of computer memory requirements. Although the method of moments has 
been repeatedly explored in the past, it is for the first time an approach has been developed to use it 
without imposing particle size distribution function. 

Also, a method was developed that allows one to combine the formalisms of the detailed 
description of the chemical nucleation of soot precursors, polycyclic aromatic hydrocarbons, and 
detailed account of coagulation and surface growth of soot particles; that is, a mathematical 
formalism which combines the linear and non-linear lumping techniques developed under the 
present program. This method was applied to analyze particle inception under shock-tube 
conditions. 


PAH FORMATION IN THE INTERSTELLAR MEDIUM 

There have been recent reports that previously unidentified the near Infra-Red emission of the 
Interstellar Medium is due to PAH. 10 It is actually estimated that PAH may be the most abundant 
organic molecules present in the Interstellar Medium. 11 In collaboration with Professor Eric 
Feigelson from the Department of Astrophysics at the Pennsylvania State University we have 
undertaken some computer simulations of PAH formation at the conditions of carbon stars using 
the mechanisms of PAH formation developed under our NASA-supported program for the 
conditions of hydrocarbon flames. The significance of such a test is the following. If our model 
can predict the observations under conditions so different from those of flames, it would further 
support the mechanism we proposed. 
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